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Abstract 

Fast rotation of compact stars (at submillisecond period) and, in particular, their stability, are sensitive to the equation 
of state (EOS) of dense matter. Recent observations of XTE J1739-285 suggest that it contains a neutron star rotating 
at 1122 Hz (Kaaret et al., 2007). At such rotational frequency the effects of rotation on star's structure are significant. 
We study the interplay of fast rotation, EOS, and gravitational mass of a submillisecond pulsar. We discuss the EOS 
dependence of spin-up to a submillisecond period, via mass accretion from a disk in a low-mass X-ray binary. 
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1. Introduction 

Neutron stars, and more generally, compact stars 
(neutron stars, hybrid hadron-quark stars, quark 
stars) are the densest stellar objects in the Uni- 
verse. Due to their compactness and strong grav- 
ity, compact stars can be very fast rotators. In the 
present paper we limit ourselves to a rigid rotation. 
Theoretical calculations show that compact stars 
could rotate at sub-millisecond p eriods (i.e., at fre- 
quency / = 1 /peri od >1000 Hz: ICook et~atlll994 ; 



Salgado et al.l ll994). 



The quest for fast rotating compact stars has 
an interesting history. The first milli second pulsar 
B1937+21, rotating at / = 641 Hz (| B acker et all 



19821 ). remained the most rapid one for 24 years. In 
January 2006, discovery of a faster pulsar J 1748- 
2446ad rotati n g at / = 716 Hz was announced 
( Hessels et all 20061) . However, sub- kHz frequen- 
cies are still too low to significantly affect the 
structure of m a ssive n eutron stars w i th M > 1M@ 
(jShapiro et all Il983t lHaensel et all . l2007t) . Actu- 



ally, pulsars B1937+21 and J1748-2446ad still ro- 



tate in a slow rotation regime, because their / is 
significantly smaller than the mass shedding (Ke- 
plerian) frequency /k- In the slow rotation regime 
rotational effects in neutron star structure, e.g., 
polar flattening, are cx (//jk) 2 -C 1, and hence not 
large. Rapid rotation regime for M > 1M© requires 
sub-millisecond pulsars with supra-kHz frequencies 
(/ > 1000 Hz). 
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Exc iting news came in December 2006. lKaaret et al 
(|2007l ) reported a discovery of oscillation frequency 
/ = 1122 Hz in an X-ra y burst from t h e X-ra y tran- 
sient, XTE J 1739-285. iKaaret et ail (|2007h wrote 
cautiously "this oscillation frequency suggests that 
XTE J1739-285 contains the fastest rotating neu- 
tron star yet found" . If confirmed, this would be the 
first detection of a sub-millisecond pulsar (discovery 
of a 0.5 ms pulsar in SN1987A remnant, announced 
in January 1989, was withdrawn one year later). 

Fast rotation of compact stars is sensitive to the 
stellar mass and to the equation of state (EOS). Hy- 
drostatic, stationary configurations of neutron stars 
rotating at given rotation frequency / form a one- 
parameter family, labeled by the central density. 
This family - a curve in the mass - equatorial radius 
plane - is limited by two instabilities. On the high 
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central density side, it is instability with respect to 
axi-symmetric perturbations, making the star col- 
lapse into a Kerr black hole. The low central den- 
sity boundary results from the mass shedding from 
the equator. In the present paper we discuss the de- 
pendence of rotation at / > 1000 Hz on the poorly 
known EOS, and derive constraints on the EOS of 
neutron stars which could result from future obser- 
vations of stably rotating sub-millisecond pulsars. 

The plan of the paper is as follows. In Sect. 1 we 
briefly describe EOSs used in our calculations. Nu- 
merical methods of solving equations of hydrostatic 
equilibrium of rigidly rotating compact stars in Gen- 
eral Relativity, and criteria for their stability, are 
presented in Sect. 2. In Sect. 3 we consider the EOS 
dependence case of rotation at / = 1122 Hz. A sys- 
tematic study of the EOS-dependcnce of fast rota- 
tion at f — 1000 — 1600 Hz is presented in Sect. 4. 
Reaching fast rotation via disk accretion in LMXBs 
and the EOS-dependence of the spin-up track is re- 
viewed in Sect. 5. Finally, Sect. 6 summarizes briefly 
main results reported in the paper. 

We are pleased to present our paper on fast ro- 
tation of neutron stars in the proceedings of the 
conference in honor of J. -P. Lasota. In 1990s Jean- 
Pierre became intrigued by a puzzling precision of 
"empirical formula" for the absolute upper bound 
on rotation frequency of compact stars, proposed 
by two of us (PH and JLZ) on the aftermath of 
i ll fated "discovery " of 0.5 ms pulsar in SN 1987A 
(|Haensel fc ZdunikL Il989h . Tw o joint papers re- 



sulted from our col laboration ( Lasota et al. . 19961 : 



Haensel et al. , 19991) , some progress has been made, 
but the basic puzzle still remains to be solved. 

2. Equations of state 

In view of a high degree of our ignorance of EOS 
of dense matter at supranuclear densities (p > 3 x 
10 14 g cm -3 ) we considered a broad set of theoreti- 
cal models. The set of ten equations of state (EOSs) 
considered in the paper is presented in Fig. [T] The 
EOSs are also listed in Table 1, where the basic in- 
formations (label of an EOS, theory of dense matter, 
reference to the original paper) are also collected. 

Two EOSs were chosen to represent a soft 
(BPAL12) and stiff (GN3) extreme cases. These two 
extreme EOSs should not be considered as "realis- 
tic" , but they are used just to "bound" the models 
from the soft and the stiff sides. 

Four EOSs are based on realistic models involving 
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Fig. 1. The equations of state in the logP - 
labels - see Table 1. 



logp plane. For 



only nucleons (FPS, BBB2, DH, APR). The next 
four EOSs are softened at high density either by 
the appearance of hyperons (GNH3, BGN1H1), or a 
phase transition (GMGS-Km, GMGS-Kp). A soft- 
ening in the latter case is clearly visible in Fig. [1] at 
pressure P ~ 10 35 dyn/cm 2 . 

Two EOSs, GMGS-Kp and GMGS-Km, describe 
nucleon matter with a first order phase transition to 
a kaon condensed state. In both cases the hadronic 
Lagrangian is the same. However, to get GMGS-Kp 
we assumed that the phase transition takes place be- 
tween two pure phases and is accompanied by a den- 
sity jump, calculated using the Maxwell construc- 
tion. The GMGS-Km EOS was obtained assuming 
that the transition occurs via a mixed state of two 
phases (Gibbs construction). A mixed state is ener- 
getically preferred when the surface tension between 
the two phases is below a certain critical value. As 
the value of the surface tension is very uncertain, we 
considered both cases. 

3. Calculation of stationary rotating 
configurations and their stability 

We computed stationary configurations of 
rigidly rotating neutron stars in the framework 
of General Relativity by solving the Einstein 
equations for stationary axi-symmetric sp a cetim e 



Bonazzola et al. . 1993 ; Gourgoulhon et al. . 1999f ). 



Numerical computations have been performed us- 



2 



Table 1 Equations of state. N - nucleons and leptons only. NH - nucleons, hyperons, and leptons. Labels of EOSs are composed 
from first letters of names of authors of EOSs of the core. In all cases but FPS, the EOS of the crust is the DH one. For FPS 
EOS its own crust model is used. 





model 


reference 


A T 1 O 
DrJ\LjLZ 


N energy density functional 


rsombaci ( lyyo ) 


FPS 


N energy density functional 


Fandnaripandc & Kavcnhall (1989) 


GN3 


N relativistic mean field 


Glendenning f 19851 


DH 


N energy density functional 


Douchin & Haensel (20011 


APR 


N variational theory a 


Akmal et al. fl998l 


BGN1H1 


NH, energy density functional 


Balberg & Gal fl997l 


GNHQm2 


NH + mixed baryon-quark state 


Glendenning ("20001 


BBB2 


N Brueckner theory 


Baldo et al. C 19971 


GNH3 


NH relativistic mean field 


Glendenning C19851 


GMGS-Km 


N + mixed nucleon-kaon condensed b 


Pons et al. (20001 


GMGS-Kp 


N + pure kaon condensed b 


Pons et al. (20001 



a A185+UIX* model of I Akmal et~al] dl998ll . 
b GM+GS model with U 1 ^ = -130 MeV. 



ing the rot star code from the LORENE library 
(http://www.lorene.obspm.fr). We calculated 
one-parameter families of stationary 2-D configura- 
tions for ten EOSs, presented in Fig.Q]and Table 1. 

Apart from fulfilling the equations of hydrostatic 
equilibrium, stationary configurations were required 
to be stable. Two instabilities were considered. 

Mass- shedding. Stability with respect to the mass 
shedding from the equator implies that at a given 
gravitational mass M the circumferential equato- 
rial radius R cq should be smaller than i? max which 
corresponds to the mass shedding (Keplerian) limit. 
The value of i? max results from the condition that 
the frequency of a test particle at circular equato- 
rial orbit of radius i? max just above the equator of 
the actual rotating star is equal to the rotational fre- 
quency of the star. This stability condition sets the 
bound on our rotating configurations from the right 
side on M — R eq plane. It fixes the largest radius 
and the smallest central density, allowed for stable 
stationary configurations. 

Axi- symmetric oscillations. Instability with respect 
to these oscillations determines the bound for most 
compact stars, with the smallest radius and the high- 
est central density (i.e., from the left side on the 
M — R eq plane). This bound is determined by the 
condition: 



For stable configurations: 



(1) 



> 



(2) 



fdhf 
V 9p c / j 

In the opposite case, {dM / dp c ) j < 0, a compact star 
star is doomed to collapse into a Kerr black hole. 

4. An example: compact stars at 1122 Hz 

In this section we present the parameters of the 
stellar configurations rotating at frequency 1122 Hz, 
a suggested rotation frequen cy of XTE J1739-28 5. 
For details and discussion see iBeiger et all (|2007ll . 

In Fig. [5] we show the M(R eq ) plots of stable 
stationary configurations rotating at / = 1122 Hz. 
We considered EOSs from Table 1. The relation be- 
tween the calculated values of M and R eq at the 
"mass shedding point" is extremely well approxi- 
mated by the formula for the orbital frequency for a 
test particle orbiting at r = R cq in the Schwarzschild 
space-time of a spherical mass M (which can be re- 
placed by a point mass M at r = 0). Let us de- 
note the orbital frequency of such a test particle 
by f^ w -(M, R cq ). The locus of points satisfying 
/ Q s r c b hw (M, R cq ) = 1122 Hz is represented by a dash 
line in Fig. [2j The points on the dash line satisfy the 
relation 
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Fig. 2. Gravitational mass, M, vs. circumferential equatorial 
radius, R e q, for neutron stars stably rotating at / = 1122 Hz, 
for ten EOSs (Fig. [TJ- Small-radius termination by filled cir- 
cle: setting-in of instability with respect to the axi-symmetric 
perturbations. Dotted segments to the left of the filled circles: 
configurations unstable with respect to those perturbations. 
Large-radius termination by an open circle: the mass-shed- 
ding instability. The mass-shedding points are very well fit- 
ted by the dashed curve _R min = 15.52 (M/1AM Q ) 1/S km. 
For further explanation see the text. 



_L ( GM \ 

2^ v^7; 



1/2 



1122 Hz 



(3) 



This formula, obtained for the Schwarzschild met- 
ric, coincides with that for Newtonian gravity for 
a point mass M. As one can see, it passes through 
(or extremely close to) the open circles denoting the 
(numerically calculated) actual mass shedding (Ke- 
plerian) configurations. This is a quite remarkable 
property, in view of rapid rotation and strong flat- 
tening of neutron star at the mass-shedding point, 
visualized in Fig. [3] 

Equation ([3]) implies a useful constraint for 
compact stars rotating at 1122 Hz: 

1/3 



15.52 



M 



1.4 M, 







km 



(4) 



5. Submillisecond pulsars 



In this section we present results for compact 
stars rotating at sub- millisecond periods (supra-kHz 
frequencies) for a broad range of frequencies, 1000 
- 1600 Hz. As in Sect. 01 calculations are performed 
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Fig. 3. Cross section in the plane passing through the ro- 
tational axis of neutron stars rotating at the mass-shedding 
limit / K = 1122 Hz (i.e., with i? cq = iJ max ), for the BGN1H1 
EOS ( z > 0) and DH EOS (z < 0). The coordinates x and 
z are defined as x = rsin0cos</>, z = rcosS, where r is ra- 
dial coordinate in the space-time metric. Dotted contour - 
crust-core boundary. 
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Fig. 4. Cross section in the plane passing through the rota- 
tional axis of neutron stars rotating at 1122 Hz at the ax- 
i-symmetric instability limit (i.e., with i?eq = R m i n ), for the 
BGN1H1 EOS ( z > 0) and DH EOS (z < 0). Notations as 
in Fig. [3] 

for ten EOSs from Table 1. Mass vs equatorial ra- 
dius relations for very fast rotating compact stars 
are presented in Fig. [5] The shape of the M(R cq ) 
curves is more and more flat as rotational frequency 
increases. For / rot = 1600 Hz the curves M(R cq ) 
are almost horizontal. At that frequency the mass 
for each EOS is quite well defined. Moreover, curves 
for different EOSs are usually well separated. Both 
features are of a great practical importance. In prin- 
ciple, they are very useful for selecting "true EOS" , 
provided one detects a very fast pulsar and is simul- 
taneously able to measure its mass. 

A dotted curve in every panel of Fig. [5] corre- 
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Fig. 5. M vs Req for stably rotating sub-millisecond com- 
pact stars. Notations as in Fig. [2] To indentify a curve cor- 
responding to a specific EOS from Table 1, one has to use 
the sequence of open circles (ordered in the same way as in 
this figure) at the mass-shedding limit in Fig. [2] 



sponds to the formula 



(5) 



used for a panel frequency (1000 Hz, 1400 Hz, 
1600 Hz). Notice that Eq. (3) is a special case of 
Eq. ([5]). Equation ([5]) works very well in very broad 
range of rotational frequencie s and EOSs (recently 
this formula has been tested bv lKrastev et al. ( 20071 ) 
for the frequency 716 Hz of PSR J1748-2446ad.) 

6. Spin up by accretion 

It is commonly believed that fast (millisecond) 
pulsars are recycled old neutron stars, spun up to 
kHz frequencies via a long-time disk accretion in 
LMXBs. Such neutron stars have a weak magnetic 
field (B < 10 10 G), which does not affect the ac- 




R eq [km] 

Fig. 6. Mass vs radius relation for an accreting neutron 
star with DH EOS. The star spins-up from / = to 
/rot = 1400 Hz. For further explanations see the text. 

cretion flow. Therefore, the accretion disk extends 
down to the innermost stable circular orbit (ISCO). 

In the present section we study some aspects of 
spin-up by accretio n from the ISCO. We use th e pre- 
scription given bv IZdunik etaD (|2002L l2005h . Spe- 
cific angular momentum per unit baryon mass of 
a matter element orbiting the neutron star at the 
ISCO, ^ ISCO , is calculated by solving exact equations 
of the orbital motion of a particle in the space-time 
pr oduced by a rotati ng neutron star ( Appendix A 



of lZdunik et alJl2002f ). 



Consider accretion of an infinitesimal amount 
of baryon mass dMe onto a rotating neutron star. As 
the star is assumed to spin up in a quasi-stationary 
manner, accretion of dMa leads to a new rigidly ro- 
tating configuration of mass Mb + dMa and angular 
momentum J + d J, with 



dJ = x t l lsco dM B ■ 



(6) 



Here, xi denotes the fraction of the angular momen- 
tum of the matter element, transferred to the star. 
The remaining fraction, 1 — xi , is assumed to be lost 
via radiation or other dissipative processes. Let us 
consider two values of xf. xi = 1 and Xi = 0.5. In 
Fig. [6]we plot the curves M(R cq ), corresponding to 
the spin-up from / = to (final) / rot = 1400 Hz. 
Calculations are performed for the DH EOS. Point 
F on this curve corresponds to the onset of instabil- 
ity with respect to axi-symmetric oscillations (con- 
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Initial mass M(f = 0) [M s ] 

Fig. 7. (Color online) The bounds on the initial mass a static 
neutron star (horizontal axis) to be spun-up up to a given 
final frequency (vertical axis). Calculations for the DH EOS. 
Labeled points (A, B, . . . ) correspond to the same points in 
Fig. \E\ The triangular region on the left (red) corresponds to 
xi = 0.5. Triangular region near the center (cyan) is obtained 
for xi = 1. The non-shaded region on the right, composed 
of a tilted triangle and a lentil-like upper determines the 
range of allowable masses of final configurations rotating at 
a given frequency. 

dition given by Eq. (Q])). Point E is the Keplerian 
configuration at frequency 1400 Hz, while point G 
corresponds to maximum mass along the curve with 
a fixed rotation frequency 1400 Hz. 

Let us consider several spin-up tracks, all termi- 
nating on the 1400 Hz mass-radius line. The curves 
starting at points A,B,C and D are the tracks of 
accreting neutron stars defined by the Eq. ( [6j , for 
xi = 1 (solid line) and xi = 0.5 (dotted line) for 
cases C,D, and A,B, respectively. In order to reach 
the frequency 1400 Hz, one has to start with a non- 
rotating neutron star located in the segment CD (if 
xi = 1), or AB (if xi = 0.5). As one can see, there 
are bounds on the initial mass of a non-rotating star, 
which can be spun-up to a given frequency / ro t via 
disk accretion. For / rot = 1400 Hz xi — 1 the allowed 
mass range for initial non-rotating star is 1.7 'Mq < 
Ah < 1.92M Q . 

In Fig [7] we plotted, for the DH EOS, the 
bounds on the initial mass of a non-rotating neu- 
tron star, for which the spin-up via accretion could 
reach a required rotation frequency. Results are 
presented for two values of xi. Shaded areas cor- 



respond to the allowed a range of initial masses of 
non-rotating star. Left "triangle" was obtained for 
xi = 0.5, while the central shaded "triangle" (with 
points C and D) - for Xi = 1. 

The limits on the actual mass of rotating NS 
which was formed by disk accretion onto an initial 
static star are given by the three curves on the right 
in Fig. [7] These curves pass by three points E,F, 
and G in Fig. [6l The curve passing by the point E 
(green) is the bound resulting from the Keplerian 
limit of a rotating star. The (magenta) line passing 
by the point F is the boundary resulting from the 
instability with respect to axi-symmetric perturba- 
tions. Finally, curve passing through the point G is 
a locus of maxima on the mass-radius curve at a 
fixed final rotation frequency: point G is obtained 
for 1400 Hz. As the frequency of the final config- 
uration increases, the mass at the Keplerian limit 
increases more rapidly, than that defined by the on- 
set of the axi-symmetric instability at the maximum 
mass (at 1400 Hz: points F and G, respectively). 
For very high frequencies the maximum mass of the 
stars rotating at a fixed frequency is given by the 
value for Keplerian limit. In the considered case of 
the DH EOS, the point G disappears at frequency ~ 
1500 Hz (point G max ). For faster rotation the curve 
M(R Qq ) is monotonous (no maximum between the 
both ends). However, the range of masses of stars 
rotating at so high frequency is very narrow. For 
/rot > 1400 Hz it is smaller than 0.1M Q (see also 
discussion of Figf5]in Sect.[5|). 

7. Discussion and conclusions 

Let us summarize main results reported in the 
present paper, and obtained using a large set of the- 
oretical EOSs. The M(R cq ) curve for / > 1400 Hz 
(rotation period < 0.7 ms) is flat. Therefore, at such 
frequencies, for any given EOS the mass of a sta- 
bly rotating compact star is quite well defined. Con- 
versely, a measured mass of a compact star rotat- 
ing at / > 1400 Hz will allow us to unveil the ac- 
tual EOS of dense matter. The "Newtonian" formula 
for the Keplerian frequency reproduces surprisingly 
well precise 2-D simulations and sets a firm upper 
limit on R eq for a given /. Finally, observation of 
/ > 1200 Hz sets stringent limits on the initial mass 
of a non-rotating star which was spun up to this fre- 
quency via accretion from a disk. 

In the present paper we limited ourselves to 
hadronic stars, built exclusively or predominantly 



G 



of hadrons. We did not discuss our results obtained 
for hypothetical compact stars built of a self-bound 
quark matter, called quark stars or strange stars. 
Some of results for quark stars with and without 



crust were briefly reported in Ref. iBeiger et al 



(2007). Generally, most of general features and rela- 
tions obtained for hadronic stars hold also for quark 
stars. However, is some cases one notices system- 
atic difference between rapidly rotating hadronic 
and quark stars. Our results will be described in a 
forthcoming publication. 
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